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1. Introduction 

Future experiments at the LHC will open a new page in the physics of hadron-hadron interac- 
tions. The new energy region promises many new physical phenomena. It is very likely that in this 
region the asymptotic properties of the hadron-hadron interaction will appear. These properties 
are connected with the first principles of the theory such as analyticity, unitarity and Lorenz in- 
variance. In the momentum transfer representation, some asymptotic bounds were obtained many 
years ago. For example, the AKM (Auberson-Kinoshita-Martin) ^ theorem leads to possible os- 
cillation at very small t in the elastic differential cross sections. Such oscillations may have already 
been discovered in the UA4/2 experiment [^, ^l. In the impact parameter representation, a more 
complicated picture appears. For example, the Froissart-Martin theorem based on the first princi- 
ples establishes a bound on the growth of total cross sections [Q] - they cannot grow faster than 
log^{s/sQ). However, the Regge behavior and the linear BFKL equation lead to a power growth of 
the total cross sections. Obviously, this conflict requires a non-linear regime at super-high energies. 
This regime may be connected with the saturation of the density of the partons in hadron and can 
be related to the saturation of the unitarity bound only in the impact parameter representation. 

The usual non-linear equations, such as the BK equation [||], are connected with perturbative 
and non-perturbative infrared effects that lead to the saturation regime at large impact parameters. 
Hence, they are not directly related to the Black Disk Limit (BDL)-saturation which in most models 
is connected with small impact parameters. Also, the models based on the color dipole represen- 
tation of the high-energy interaction [^] suppress, for the most part, the soft hadron interaction. 
This is the reason why these models describes the hard interactions well. On the other hand, the 
BDL suppresses in most cases the hard component. It leads to an increase in the role of the soft 
component at super-high energies. 

Unitarity of the scattering matrix, 55+ = 1, is connected with the properties of the scat- 
tering amplitude in the impact-parameter representation as the latter is equivalent at high energy 
to a decomposition into partial-wave amplitudes. At high energies the scattering amplitude in 
the impact-parameter representation can saturate the unitarity bound of the overlapping function 
G{s,b) below some impact parameter bi. One can impose the unitarity condition through various 
schemes. Two of them are based on an analytic implementation of unitarity: the first is the standard 
eikonal representation 

poo 

T{s,t) ~ / / bdbJoibq) [1 - txp{-x{s,b)], (1.1) 

with t = —q^. The second corresponds to the ?7— matrix approach [^: 

G{s,b)=U{s,b)/[\ -iU{s,b)]. (1.2) 

where U {s,b) = ix{s,b) = Thom{^,b). In [||] it was shown that such a solution can lead to a non- 
standard behaviour of the ratio a^/ / Gtot —>■ 1 , as 5 ^ oo. 

These schemes lead to a decrease in the growth of Gtot, but they predict very different energy 
regions for the onset of saturation. We show also that such a decrease heavily depends on 
the profile in impact parameter dependence of the Born amplitude. In particular, if one takes the 
eikonal phase in the factorized form, one usually supposes that despite the fact that the energy 
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dependence of h{s) can be a power of h{s) ~ s^, the total cross section will satisfy the Froissart 
bound Ofot < ct log^{s). We find in fact that the energy dependence of the imaginary part of the 
amplitude and hence of the total cross section depends on the form of f{b), i.e., on the s and t 
dependence of the slope of the elastic scattering amplitude. For example, a power dependence on 
the impact parameter will still violate the Froissart-Martin bound after eikonalization. 

1.1 Saturation and non-linear equations 

Different approaches of the saturation of the scattering amplitude try to connect it with the 
non-linear saturation processes which have been considered in a perturbative QCD context [10, ll|]. 



Such processes lead to an infinite set of coupled evolution equations in energy for the correlation 
functions of the multiple Wilson lines. In the approximation where the correlation functions for 
more than two Wilson Unes factorize, the problem reduces to the non-Unear Balitsky-Kovchegov 
(BK) equation 

It is unclear how to extend these results to the non-perturbative region, but one will probably 
obtain a similar equation. In fact we found simple differential equations that reproduce either the 
U -matrix or the eikonal representation. 

Let us consider the saturation equation that has been used for many years in the various fields 
of physics - that is the logistic equation: 

dN/dy = dN {l-XN); (1.3) 

This equation has a standard solution in the form 

= f{b) e'' ^ y/[\ + A f{b) e^'^y]. (1.4) 

If, inspired by the BK results, we take d = a{0) — I with a(0) the intercept of the Pomeron and use 



the evolution variable y = logs, with A = 1, Eq. (L4) coincides with the standard ?7-matrix rep- 
resentation of the scattering amplitude obtained in the ?7-matrix formalism for 3U {s,b) = s^f{b). 



Hence, Eq. (1.4) extends the unitarized form of the scattering amplitude and is determined by 



the non-linear equation ( |1.3[ ). We can see that the asymptotic regime (y = log (5) 00) leads to a 



bound 1/A for the amplitude, and gives the standard BDL bound for A = 1. We can use the same 



procedure [12] to obtain other unitarized forms of the scattering amplitude. 



For example, from the non-linear equation 

dN/dy = d{\ -XN) log(l - AA^), (1.5) 

we obtain a solution 

N=U\ -e-''lf^\ (1.6) 
A 



The equation reproduces the extended eikonal (quasieikonal) form [[l3|]. With A = 1, Eq. (1.6) 



coincides with the standard eikonal form of the scattering amplitude in the impact parameter rep- 
resentation. 

We show in Fig. 1 and Fig. 2 the dependence of these non-linear equations on the size of the 
intercept of the Pomeron and on the energy S = ^/s in the case A = 1. We see that the qualitative 
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Figure 1: dN/dy for Eq. (1.2) - left - and for Eq. (1.5) - right -. 



behavior of both equations is very similar. The derivative in Eq. ( [T^ ) has a larger maximum at low 
s than that in Eq. dOj ). However, at large s, the situation in the opposite: the derivative in Eq. 
exceeds that of Eq. ( [L6| ). 

We have shown that the most usual unitarization schemes could be recast into differential 
equations which are reminiscent of saturation equations [^. Such an approach can be used to 
build new unitarization schemes and may also shed light on the physical processes underlying the 
saturation regime. 




Figure 2: dN/dy for Eq. (1.2) and for the Eq. ( 1.5) 
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As an example of other unitarization schemes, the equation 

dNbare/dt[l-N^] (1.7) 



has the unitarized solution 



1.2 Unitarization 



dy 



G{s,b) = tanh[A^fo^,,]. (1.8) 



Let us now compare the behavior of the scattering amplitudes in different unitarization schemes. 
They are given by 

poo 

F{s,t) = i bdbG{s,b) Jo{qb), (1.9) 
Jo 

where G{s,b) is the overlapping function which depends on the scheme and satisfies the BDL 
bound G{s,b) < 1. G{s,b) is a function of s, b, Nbare< and 

/•OO 

Nbare{s,b) =i q dq FBom{s,q) Joiq b). (1.10) 

JO 

If the Born scattering amplitude is in the Regge form 

FBorn{sj) = i h exp{R^ /I t) (1.11) 
the corresponding form in the impact parameter representation is given by 

FBom{s,b) = ih ^ exp{-b^/{AR^)). (1.12) 

This Bom amplitude will be the input of the different unitarization schemes. 

Following the Van Hove normalization 0], the corresponding total cross section is 

noo 

o,o,{s) =An b[G{s,b)]db. (1.13) 

JO 

In Fig. 3, we show the total cross section for different unitarization schemes. For all the 
schemes, the Born amplitude was taken in the simplest Regge form ( 1.11 ) with the parameters 
corrsponding to the exchange of a hard Pomeron: 

h = 0.1; R = 3GeV"^ d = a(0) - 1 = 0.4 

Here we want to compare the schemes directly. It is clear that if we make a fit to experimental data, 
the constant h can be different for different unitarization schemes. 

We can see from Fig. 3 that at small energies {^/s = 10 GeV) the total cross sections are prac- 
tically the same for all unitarization schemes. At this energy the unitarized amplitude practically 
reduces to the Born amplitude. When the energy grows, the schemes differ most in an intermediate 
regime which corresponds to the energy region of approximately 30 GeV < < 200 GeV. 
This regime corresponds to the growth of the overlapping function which is different in different 
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Figure 3: Total cross sections for different unitarization schemes. The hard, long-dashed and short-dashed 
lines correspond to the eikonal, the f/-matrix and the tanh forms of unitarization. 



unitarization schemes. However, when we come to the asymptotic regime (with saturation of the 
BDL bound) the difference ceases and the growth is practically the same in a wide energy region. 

Hence, the asymptotics of all unitarization schemes are practically the same. Such a behavior 
leads to two methods of analysing experimental data. One is connected with a different contribution 
of the secondary Reggeons to the cross section. In this case, the asymptotics of the total cross 
sections can be exactly the same for different unitarization procedures. The second solution is 
obtained if we take the same contributions of the secondary Reggeons. The asymptotics of Otot 
will then have the same energy dependence but some quantitative difference between the different 
unitarization schemes. In this case, the comparison of the predictions for the total cross sections at 
LHC energies from the different fitting procedures with different unitarization schemes will point 
out the true unitarization scheme and the true form of the non-linear equation. 



2. Conclusion 



In the presence of a hard Pomeron [15], the saturation effects can change the behavior of some 
features of the cross sections already at LHC energies. Hence, the problem of unitarization of the 
scattering amplitude will be very important for the analysis of experimental data. 

Non-linear effects which work in the whole energy region supply an acceptable growth of 
the total cross sections. Saturation leads to a relative growth of the contribution of peripheral 
interactions. The most usual unitarization schemes can be recast into differential equations which 
are reminiscent of saturation equations. We can see that the different forms of the non-linear 
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equations lead to the different behavior of the scattering ampUtudes at intermediate energies and 
to the same behavior at asymptotic energies. Further exploration of the unitarization schemes and 
of their connection with the future experimental data of LHC may shed some light on the physical 
processes underlying the saturation regime at the partonic level. 
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